Decoherence and Entanglement Dynamics of Coupled Qubits 
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We study the entanglement dynamics and relaxation properties of a system of two interacting qubits in the 
two cases (I) two independent bosonic baths and (II) one common bath, at temperature T. The entanglement 
dynamics is studied in terms of the concurrence C (i) between the two spins and of the von Neumann entropy 
S{t) with respect to the bath, as a function of time. We prove that the system does thermalize. In the case (II) of 
a single bath, the existence of a decoherence-free (DFS) subspace makes entanglement dynamics very rich. We 
show that when the system is initially in a state with a component in the DFS the relaxation time is surprisingly 
long, showing the existence of semi-decoherence free subspaces. The equilibrium state in this case is not the 
Gibbs state. The entanglement dynamics for the single bath case is also studied as a function of temperature, 
coupling strength with the environment and strength of tunneling coupling. The case of the mixed state is finally 
shown and discussed. 
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I. INTRODUCTION 

The principle of superposition is the most important fea- 
ture of quantum mechanics. It gives rise to interference and 
to quantum entanglement. Among its endless implications, is 
the possibility of processing information in a way that goes 
beyond the classical scheme of the Turing machine, with the 
fundamental consequence that information is physical. Com- 
puter science is thus not a branch of pure mathematics. More- 
over, quantum computation promises to be able to solve cer- 
tain computational problems with an exponential speed-up 
(for a review on quantum information, see fv\). 

A quantum system will generically lose coherence when in- 
teracting with another quantum system. In this case, its evolu- 
tion will not be unitary, and there will be quantum correlations 
with the environment. The phenomenon of decoherence has 
been advocated as the solution of the measurement problem in 
uantum mechanics, and the appearance of the classical world 
[bo . For quantum information, preserving coherence is the 
most important and demanding problem in order to build a 
functioning quantum computer [4J. The loss of quantum be- 
havior in a system can be measured by the loss of quantum 
correlations within the system, and the increase of quantum 
correlations with the environment, or, in other words, by the 
entanglement within the system, and with the environment. 
Unfortunately, there is not a general way to study entangle- 
ment between three parties. A relevant exception is consti- 
tuted by the case of two qubits because the concurrence C be- 
tween the two qubits is a valid measure of their entanglement 
even in a mixed state |5]. Therefore, in the case of a compos- 
ite quantum system constituted by two qubits and an environ- 
ment, the entanglement between the qubits can be character- 
ized by their concurrence, while the von Neumann entropy S 
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can measure the entanglement between the system of the two 
qubits and the environment. 

In this paper we study a quantum system made of two in- 
teracting qubits coupled to a bosonic bath. We consider the 
following two cases: (I) for each spin one independent bath 
and (II) both spins coupled to the same bath. 

The non equilibrium dynamics of the reduced density ma- 
trix p(t) of the two qubits is studied within the Bloch-Redfield 
approach [6]. We assume that the environment is in thermal 
equilibrium at temperature T and the system-bath coupling 
is weak. With this assumption, one obtains a set of coupled 
integro-differential equations for the elements of the reduced 
density matrix. In the Markov approximation, the equations of 
motion for p{t) take then the form of simple linear differential 
equations that can be solved analytically. Although quantum 
systems interacting with a bath are generically believed to re- 
lax, there are only a few examples where one can carry out all 
the calculations and prove directly how relaxation is achieved. 
In the case (I) we prove that the system always relaxes to the 
Gibbs state. The case (II) is more rich due to the existence of 
a decoherence-free subspace (DFS) J13t]. The novelty of our 
approach is also in taking in explicit consideration what hap- 
pens if the system is initially in a state which has a component 
in the DFS, which has very important consequences. 

The main object of this work is the study of the quantities 
C{t) and S{t) as a function of time and temperature T for 
both the cases (I-II). Since we are interested in the formation 
of quantum correlations, we put emphasis in the case T = 
where the entropy S has the meaning of entanglement with the 
bath degrees of freedom. The system of two spins in bosonic 
bath has been studied in several papers 1171 I8ll9l [l0l [Till . 

As mentioned above, the existence of the DFS reveals an 
important novel effect. There is an "interference" effect that 
decreases dramatically the decoherence rate when the initial 
state is in a coherent superposition between the DFS and its 
complement. We also see that the system relaxes to a state that 
is not the Gibbs state because the amplitude of the singlet state 



must stay constant. Moreover, we study entanglement dynam- 
ics and relaxation behavior for the system (II) as a function of 
temperature T and coupling strength. Finally, we study the 
case of the mixed state, to show that the dramatic increase of 
relaxation time discussed above is genuinely a quantum effect. 



n. THE MODEL 

We consider a simple generalization of the spin-boson 
problem 11411 . where two qubits interact with each other via 
an Ising type coupling and are also coupled to a bosonic envi- 
ronment. The system Hamiltonian is 



Hs 



A A 



2 2 2 

We take h = ks ~ I- Here ct's and r's are the Pauli matri- 
ces on the first and the second spin respectively. For simplicity 
bias terms are absent and the tunneling coupling A is the same 
for both spins. The system Hamiltonian is trivially diagonal- 
ized (see appendix) and we will denote by {|i?i)} the basis 
of its eigenstates. In case (I), the baths are modeled by the 
Hamiltonian 



i? 



(2) 



The interaction Hamiltonian is 



H^ 
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1 



fT.^C, 



■«a(^i,l+^".l)+9'^^Xl^"^2(&i,2 + ^a,2) 



(3) 

Each spin is subject to one fluctuating force and the two forces 
are uncorrected. The coefficients Ca,i give the strength of 
the coupling of the system to each harmonic oscillator of the 
bath, which we assume to be weak. For a gaussian model 
all the properties of the baths are described by the spectral 
density Ji{uj) — tt^^ c^ ^6{u!a^i — lo) which, in the case of 
ohmic baths that we consider here, takes the form Jiiuj) = 
2TiKi LO exp(— w/wc). Here lUc is the cut-off frequency which 
is assumed to be the largest energy scale in the problem. We 
study the case Ki = K2 = K/27r, with the two baths at the 
same temperature T. In the case (II) the same fluctuating force 
acts on both spins so that Hb has only one set of harmonic 
oscillators. The coupling of the spins to the reservoir in this 
second case reads 
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III. BLOCH-REDFIELD APPROACH 

In this section we introduce the equations which allow us 
to study the time evolution of the reduced density matrix of 
the two spins. Within a Markov approximation, the matrix el- 
ements of the reduced density matrix in the eigenvectors basis 
{|£'i)} obey the following generalized master equation ||6|]: 
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is the Redfield tensor defined by R. 
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The rates F^ are given in the appendix. The indexes 
m,m'n,n' run from 1 to 4. Ei,...,E4 are the eigenvalues 
of Hs and cj„i„ = Em — En ■ 

Equation ^ can be simplified by keeping only secular 
terms, i.e. those terms such that the argument of the expo- 
nential is zero. The equation for the diagonal elements (popu- 
lations) of the density matrix is 

Pmm{t) = E Pnn{t)Wmn ~ Pm7n{t) ^ W„m , (6) 

with Wmn = ^tm-mn + Krnmn = 2Re r+„„„„ For the off 



(1) diagonal elements (coherences) we have 



Pmnit) = { — i'^mn — lmn)Pmn{t) + ^ Qmnkl Pkl{t), (7) 
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With 7„,'m — l^k\^ m'kkrn' ^ ^ mkkm) ^mmm'm 

^mmm'Tn' ^"^ ^mnkl = ^ Inrnk + ^ Inmk ^^ ^""i = ^kl and 

zero otherwise. Due to the symmetries of the model the coef- 
ficients ^mmm'm' ^^ identically zero. To simplify the calcu- 
lation we disregard the Lamb shift to the eigenfrequencies due 
to the imaginary part of the F's. Within this approximation we 
can write 



1 
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{Wkrn' + Wkm) . 



(8) 



We solve Eqs. (|6|l and O by Laplace transform and discuss 
separately the cases (a-b). 

Direct calculation shows that all the coefficients which ap- 
pear in Eqs. (|6]l and O are not independent. Before we pro- 
ceed to brute force solution of the differential equations it is 
convenient to employ all the symmetries of the problem. From 
the spectrum of Hs we obtain that ui^i = LO42 and W43 = W21. 
Moreover Wm„ = eyi\)[~ j3ujrnn\Wnm- For both cases (I) and 
(II) the solution of Eq.(|7]i is straightforward once we notice 
that for this model all the 8's are zero. This simplification 
does not apply when the two spins have a different A's and 
one hads to take care of all term appearing in EQ^(|7]). We 
obtain 



Pi] [i) = Pij (0) exp[- {-fij + iujij )t\ i ^ j. 



(9) 



In the case (I) direct calculation shows that VF43 = W21 and 
W42 — W31 and W41 = W32 — 0. In the case (II) the system 
is assumed to interact with a single ensemble of harmonic os- 
cillators. It is well known that in this case a DFS exists, i.e., a 
subspace of the non interacting system which is fully decou- 
pled for the environment and where the time evolution is uni- 
tary fiy\. For this simple model the DFS is one-dimensional 
and coincides with the singlet eigenstate |3). It follows that 
Wis, W^23 5 W43 are all identically zero. Moreover also 1^14 
is zero as the straightforward calculation shows. Also in this 
case we are reduced with only two parameters VF42 and W21 
which in this case are not equal as in the double bath case. Of 
course the population P33 of the singlet state remains constant 
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FIG. 1: (Color online) Entanglement dynamics for the system with 
A — l,v = 0.7 for the case (I) of two independent baths at low 
temperature P = 10. All the graphs show the Von Neumann entropy 
of the system of two spins and Concurrence between the two spins 
as a function of the time t. The initial states are |^a) — l/%/2(| Ti 
) + I iT)) G n± and 1*6) = I TT) G Hx, |*c) = i3) £ H^ 
\^d} = I Ti) G Ti., respectively. The two spins are initially disen- 
tangled with the environment. At t = oo, the system always thermal- 
izes in the ground state and therefore is a pure state: at short times the 
system gets quickly entangled with the environment due to fast deco- 
herence, then it dissipates to reach the ground state. The relaxation 
time is of the order of t ~ 400. Notice the different behavior of the 
concurrence in the different cases. In particular, there are no oscilla- 
tions for the case of the initial singlet state. When the system is very 
entangled with the environment, the concurrence dynamics changes 
qualitatively, before dissipation intervenes to damp its oscillations. 



FIG. 2: (Color online) Entanglement dynamics for the system with 
A — l,v — 0.7 for the case (II) of a single bath at low tem- 
perature l3 = 10. All the graphs show the Von Neumann en- 
tropy of the system of two spins (blue) and Concurrence between 
the two spins (red) as a function of the time t. The initial state is 
i*a) = l/\/2(| Ti) + I it)) G H±. The two spins are initially 
disentangled with the bath. At f — oo, the system thermalizes in 
the ground state and thus is disentangled again. At short times the 
system gets quickly entangled with the environment due to fast de- 
coherence, then it dissipates to reach the ground state. Unlike the 
double bath case, there is an important relationship between entan- 
glement with the bath and concurrence. The concurrence decreases, 
but when the system is very entangled with the bath, the structure of 
the oscillations changes qualitatively before the dissipation becomes 
relevant. The relaxation time is of the order of f ~ 400. 



in time. In both cases (a-b), the analytic solution of Eq.© 
is then obtained by taking the Laplace transform with respect 
to the time and by solving a set of algebraic equations. The 
explicit expressions are reported in the appendix. We can 
write the full solution of Eqs.(|6l|7|i in a superoperator form: 
p{t) = £tp{0). 



IV. ENTANGLEMENT DYNAMICS AND THE DFS 
INTERFERENCE 

Having found the solution p{t) for the time evolution of the 
density matrix of the system, we can proceed to the study of 
the entanglement dynamics. Let us call Ti" the subspace of 
the singlet state |3) and 7^-'- = span{|l), |2), |4)} the sub- 
space orthogonal to it. We will initialize the system in the 
pure state p(0) = |^) (^| corresponding to a generic superpo- 
sition between the two subspaces; j^P) = A\(j)±) + B\3), with 
|(/)_l) € H^. Having obtained the solution p{t) for the dynam- 
ics of the reduced system, we can study the von Neumann en- 
tropy S{p{t)) = — Tr(plogp) and the two-spins concurrence 
C{t) defined in 15|] 
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FIG. 3: (Color online) Here we show the entanglement dynamics 
for the same value of parameters of the previous figure but for the 
initial state \'^t) = | |T) G Ti.±. Since also this state does not have 
any component on the DFS Tis, the qualitative features of 5(f), C(t) 
are the same than for the initial state \'^a) though the details of the 
evolution are different. The relaxation time is here t ~ 400 as well. 



C{p{t)) = max(0, VA 



A4), (10) 
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FIG. 4: (Color online) Here we show the entanglement dynamics for 
the same system of Fig. 2 but with initial state \^d) = | Ti)- The 
component in H," gives a dramatically longer relaxation time. The 
Von Neumann entropy will not equilibrate to zero so that there will 
be a residual significative entanglement with the bath. The concur- 
rence oscillations are damping very slowly. In the inset, we show the 
behavior of C(i) and 5(f) at a very large f. The relaxation time is of 
the order of i ~ 4 x 10^. Moreover, the system-bath entanglement 
(Von Neumann entropy) reaches a minimum before slowing rising to 
the final value at the equilibrium. 



as a function of time. Here A^'s are the eigenvalues of picfy ® 
Cy)p* {o'y®<^y)- We have performed a study in the space of the 
parameters u, A, k, /3 = 1/T and different initial states p(0). 
The strength of the coupling n only changes the time scale 
of the evolution. The parameters v, A both set the time scale 
and a temperature scale. Moreover, for the case of the single 
bath they have an effect on the concurrence at the equilibrium. 
The effect of the temperature is that of making the system 
dissipating faster, and mixing the system. Here, we want to 
focus on the important case of zero temperature because at 
T = 0, the entropy S{t) has the meaning of measuring the 
entanglement with the bath degrees of freedom. At T > 0, 
it is not possible to distinguish the quantum correlations with 
the bath and the mixing due to the finite temperature. Initially 
the system is prepared in a pure state, and therefore S'(O) = 0. 
Then, the dynamics given by the interaction Hamiltonian will 
entangle the system with the environment. Let us examine 
now the two cases. 

(I) Two independent baths. — In Fig[T]are plotted C{t) and 
Sit) for the four initial states |*„) = 1/^2(1 Ti) + I it 
)), Y^b) = I TT) in the irreducible subspace Tij^, the singlet 
state \^,) = |3> = 1/V2(| Ti> - I it)), and 1*^) - | ti). 
The choice of parameters is w — 0.7, A — 1,k, — 0.01 and 
j3 — IQ that is a very low temperature for the system. The con- 
cuiTence C{t) goes to zero as the system decoheres. Since the 
time of decoherence is much faster than dissipation, most of 
the concurrence is lost at earlier times, then it degrades more 
slowly. If the system is initialized in Ti^, it will be oscillat- 
ing with damping oscillations. If it is initialized in the singlet 



state, it just leaks towards the other states (in particular, the 
ground state) without oscillations because |3) is an eigenstate 
of Hs (FiglT};). In all cases the concurrence revives after hav- 
ing gone to zero, in the point where the entanglement with the 
environment is maximum, a sign of monogamy of entangle- 
ment. An important effect is that, when the system gets more 
entangled with the baths, the concurrence dynamics changes 
qualitatively and the some oscillations increase their ampli- 
tude (see, in particular, FiglTJ)). 



We can prove that the system relaxes taking the limit 
Peq = liint^oo p'^''*)(i), and obtain, for a generic /3, 



H^oo P'"' (1^)7 aiiu ouiaiii, lor a generic 
p-' = Z-ldiag(e-'3'^^e-^^^e-'^^^e-'3^*) where Z == 
J2i=i e"''^* is the partition function of the system. The sys- 
tem thermalizes in the Gibbs state. This is what everyone 
would expect, a system is supposed to thermalize in the Gibbs 
state. It is though always very difficult to prove relaxation to 
equilibrium in concrete examples and this is one of our results. 
At T = the Gibbs state is the ground state and therefore the 
asymptotic value for S is zero. At T > 0, the von Neumann 
entropy will first increase to a high value due to the fast deco- 
herence, and then, when dissipation kicks in, slowly decrease 
to its asymptotic value given by the entropy of the Gibbs state. 
Notice that the relaxation time is of the order of t ^ 400, in- 
dependently of the initial state. 

(II) Single bath. — The behavior for the system in a single 
bath is completely different. In this case, the singlet subspace 
W is a DPS | 131 . If the system is prepared in the singlet state 
p{0) — |3)(3|, it will stay there forever, so the case yl = is 
trivial. If we prepare the state in the subspace T-C±, there can 
be no effect due to the DPS. Nevertheless, the entanglement 
dynamics is very interesting. In Pig|2^, the concurrence C{t) 
and the von Neumann entropy S{t) are plotted for the initial 
states |*a) = 1/V2(| ti) + I it)) e Wx- At the beginning, 
the system is not entangled with the bath, and the concurrence 
starts decreasing, due to fast decoherence, in a fashion very 
similar to the case of the double bath. Then notice, that also 
here, when the system is very entangled with the bath, the 
concuiTence dynamics changes again qualitatively. When the 
dissipation becomes relevant, the oscillations damp and the 
system relaxes to the Gibbs state. At the plotted temperature 
j3 — 10, this is practically the ground state so that the entan- 
glement with the bath is zero. The equilibration time is of the 
same order of magnitude than the case (I) with the two in- 
dependent baths. In PigEJ) the initial state is | '!'(;) = | ti), 
which is the equal superposition between I'i'a) and the singlet 
state l^c) = |3). Now things change dramatically. The con- 
currence dynamics is qualitatively completely different, and 
it damps in an extremely slower way. In Pig|4]is shown the 
behavior at large times: the equilibration time is of the or- 
der of f ^ 4 X 10^, three orders of magnitude more than 
the usual. How is it possible? Why the oscillations of the 
system are damped in such a slow way? Pven if "half of 
the system is in the DPS, it could still be that the decoher- 
ence time will depend only on the part that is in T-i±. This 
is exactly the case if we prepare the initial state in the mixed 
state p(0) = l/2(|*i)(*i| + |3)(3|). The evolution equations 
p{t) — £tp are linear and we would expect the same type of 



oscillations than for the first term only, since the second one 
is constant in time. But that is not the case. The state is pre- 
pared in a coherent superposition of |\l/i) and the singlet state 
1 3), so this means that the off-diagonal terms in p(0) show up 
in p{t) — £tp{0) in a way that makes the decoherence much 
slower. We say it is an interference effect because it discrim- 
inates between a coherent superposition and a classical mix- 
ture. As we shall see in section IVll there is no such effect if 
the initial state p{Q) is a classical mixture of states in Ti^ and 

In order to understand the long decay time when the system 
is initialized in a state with a component in the DFS Ti", let us 
look at the decay of pi3 and compare it to the one of, say, pi2. 
The oscillation behaviors are comparable but the exponential 
decays are very different. These decays are governed by the 
rates 713 and 712, as one can see from Eq.©. Using Eq.® 
and the fact the Wi^ is zero we readily find that at low temper- 
ature 713 ex coth(/3a;2i/2) — l while 712 ex coth(/3a;2i/2) + l. 
With our choice of the parameters we have 712/713 ~ 10^ as 
illustrated in the plots. 

Not only the system thermalizes with a longer time scale, 
but it does not relax to the Gibbs state. After all, we expect the 
population of the singlet state to be a constant of the motion. 
Taking in p("°"°)(i) the Hmit for t ^ cx), we find 



(mono) 

Peq 



where Z3 



diag 



\Al 



-pEi 



\A? 



-PE2 



\B\ 



\A\' 



-PE4 



(11) 



3 = X^Ms "S ''^' is the partition function over the 
irreducible subspace 7i^. Because of the constant term 

P:a° ~ I^P' the entanglement with the bath can be non 
vanishing even at the equilibrium at zero temperature. For in- 
stance, at T = the equilibrium value for the Von Neumann 
entropy S is given by 5" = \A\^ log \A\'^ + (1 - \A\^) log(l - 
\A\'^) and is obviously maximized by \A\'^ = 1/2 for which 
we have S* = 1, as it was also argued in \\Q\. The concur- 
rence C at zero temperature still depends also on w, A. For 
the case study ofv = 0.7, A = 1 the equilibrium concurrence 



at T = is Ce, 



0.33. 



V. STUDY IN TEMPERATURE AND COUPLING 
STRENGTHS 

In this section we study the behavior of the entanglement 
dynamics in the system (II) as a function of the coupling 
strength with the environment k, the temperature (3, and the 
pai-ameter A. The initial state is |*q) = l/\/2(| ti) + | it))- 

The results of the study in n are plotted in Fig.©, for a low 
temperature /3 = 10. We expect that a greater coupling with 
the environment will not make qualitative changes, as long as 
the hypotheses of weak coupling are still satisfied. We expect 
that the time scale of the system will shrink for larger cou- 
plings. Fig.® confirms this physical insight in the behavior 
of both C{t) and S{t) at every time scale. 

The second study is the behavior of the system in temper- 
ature. In Fig.® we have plotted the time evolution of S{t) 
and C{t) for different temperatures /3. At high temperature 
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FIG. 5: (Color online) Von Neumann entropy S{t) (a) and Concur- 
rence C(t) (b) as a function of time for different values of the cou- 
pling with the environment k = 0.1, 0.2, 0.3, 0.4, for the initial state 
is i*a) = 1/V2(| Ti) + I it))- The temperature is /3 = 10 and 
V — 0.7, A = 1. A stronger coupling has the effect of shrinking the 
time scale for the entanglement dynamics without any other qualita- 
tive changes. 



/3 = .1, the system decoheres very rapidly and the entangle- 
ment dynamics is trivial: it entangles and mixes with the en- 
vironment and the two spins disentangle from each other At 
a medium temperature /3 = 1, the process of entanglement 
and dissipation towards the environment is smoother, and the 
concurrence dynamics is less trivial, eventually though, the 
two spins disentangle from each other. At low temperatures 
P — 5, 20, the system shows the most interesting behav- 
ior. Now the Von Neumann entropy can be interpreted as 
just the entanglement with the environment. Its rising and 
decay marks the two phases of decoherence and dissipation. 
Initially, the system rapidly decoheres by entangling with the 
environment. The time scale of decoherence is much smaller 
than the one of dissipation. Then, dissipation intervenes, and 
at very low temperature the system must fall into the ground 
state, which is a pure state and thus disentangled with the en- 
vironment. 

The behavior of ConcuiTence C{t) is the most interesting 
as a function of temperature. At high temperatures, the con- 
currence is rapidly damped down because of the entanglement 
with the environment. The entanglement is monogamous so 
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FIG. 6: (Color online) Study in temperature of the entanglement dy- 
namics for the case (II) of a single common bath. Here v — 0.7, A = 
1,K = 0.01. The initial state is |*„) = 1/^2(1 U) + I it))- In 
the graph (a) is plotted the Von Neumann entropy S{t) as a function 
of time for different values of the temperature /3 — 20, 5, 1, 0.1. At 
low temperatures the state approaches the ground state as equilibrium 
state. The non monotonic behavior of S{t) shows the different time 
scales of decoherence and dissipation. In the graph (b) we show C(t) 
for the same temperatures. At low temperatures the phenomenon of 
entanglement revival occurs. 



as the system entangles with the environment, the mutual en- 
tanglement between the spins decreases. Once zero, even 
though the system gets more and more mixed, the Concur- 
rence remains zero. At high temperature there is no quantum 
correlation left in the system. At low temperatures instead, the 
ConcuiTence revives and its behavior is non monotonic. The 
ConcuiTence equilibrates to a non-zero value at 

with damped oscillations. The amplitude of the oscillations 
changes pattern when the system is very entangled with the 
environment. 



is at the extremely low temperature (3 — 20. Let us define 
the density matrices pa = |*a)(*a| = (| Ti> + I iT>)((Ti 
I + (it |)/2, Pb = \^b){^b\ - I TT)(TT I, Pc = |*c>(*c| = 
|3)(3| = (I Ti) - I man l - (it l)/2. we want to show 
that if we prepare the system in the initial classical mixture 

Pmix — (Pa + Pc)/'2', there is no trace of the behavior ob- 
tained when we have a coherent superposition of a state in the 
DFS Hs with one in H^. We will compare the entanglement 
dynamics with the initial preparation of the classical mixture 
of two states in Ti^, namely p^J-^ = {pa + Pb)/'^- The plot of 
Figl?] shows the results. Let us first look at the evolution for 
the entanglement S{t) with the bath. The state pU,^^ gets very 
entangled with the bath during the fast decoherence period. 
Then when dissipation becomes important, the entanglement 
decreases. Nevertheless, the equilibrium state cannot be the 
Gibbs state because the population in the DFS Ti" is constant. 
Therefore the final state is not the ground state for the sys- 
tem and some mixture is present. This mixture does not mean 
that there is residual entanglement with the bath even though 

(2) 

we are at extremely low temperature. The state p)Ji^^ instead, 
dissipates towards the ground state because there is no initial 
population in the initial state and it is, therefore, a pure state 
at the equilibrium, disentangled with respect to the bath. In 
both cases, the relaxation times are comparable t ~ 400 like 
in the case of Figs|2][3] The concurrence C{t) shows a similar 
pattern. In both cases we have revival of the concurrence after 
it hits zero and similar equilibration values. Again, the graphs 
show that the relaxation time is about t ~ 400. 

This study shows that the dramatic increase in the relax- 
ation time shown in Fig|4] is due to the quantum superpo- 
sition of state in Hs and Ti.-^. It is a purely quantum ef- 
fect. This means that if we consider for instance the subspace 
j^serrn ^ span{|3), |*a)}, although it is not a DFS, the re- 
laxation time for states in this subspace is much larger than 
that of states its orthogonal complement, and this is because 
j^serm contains a DFS. We call H'"'"' a semi decoherence- 
free subspace. The existence of such subspaces is important in 
quantum computation because one can protect quantum mem- 
ory and quantum information processing even in absence of a 
real DFS, which is often the case in presence of perturbations. 
Moreover, it enlarges the dimension of the subspace in which 
the information is protected. In the case studied for instance, 
the DFS is trivial because is one-dimensional, and no infor- 
mation (classical or quantum) can be encoded. Nevertheless, 
the subspace 7^"^'"* is two dimensional and one can encode 
a qubit in it. So an aiTay of pairs of spins 1/2 could consti- 
tute a good quantum register, of course for the model of noise 
presented here. 



VII. CONCLUSIONS AND OUTLOOK 



VI. THE MIXED STATE CASE 

In this section, we study the behavior of the system pre- 
pared in a initial classical mixed state. The common bath 



In this article, we have thoroughly studied the system of two 
interacting qubits in a bosonic environment, in both the cases 
where each spin is interacting with its own bath, and where 
they share a common bath. We solved the master equation for 
the reduced system of the two qubits and studied the entangle- 
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FIG. 7: (Color online) S{t) and C{t), for low temperature (3 = 20. 
The values for the other parameters are v = 0.7, A = 1, k = 0.01. 
The initial states are (a) the classical mixture p^^^ of two states be- 
longing to Hs and 7i^. In (b) we see the entanglement dynamics 
for the state p ^^ which is a classical mixture of two states belong- 
ing both to Ti . The relaxation times are similar. Both states show 
revival in the Concurrence. The final entanglement with the bath is 
different in the two cases because pmix cannot be the ground state 
since in the mixture the component on Tis is constant while Pmia; 
goes to the ground state. 



in the Gibbs state. Upon the completion of this work, we no- 
ticed the study of 1 15], showing, for a different model of two 
qubits in two baths, similar results for the concurrence dynam- 
ics. 
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APPENDIX 

1. Eigenvectors and eigenvalues of Hs 



We take as basis the eigenvectors of azTz'. 
— h), I ). In this basis Hs simply reads 



Hs 



We have the following eigenenergies 
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A V 



El = --y^v^+AA^ E2 



The corresponding eigenvectors are 



|1) = [r+.s+,s+,r+] 



ment dynamics by means of ConcuiTence and von Neumann 
entropy. In the case of the single bath, the entanglement dy- 
namics is very rich. The concurrence oscillates violently in 
a time window when the system is very entangled with the 
environment but the dissipative effects are not yet important. 
If the state is initialized in a coherent superposition with a 
component in a decoherence-free subspace, there is an inter- 
ference effect that changes the decoherence rate and the en- 
tanglement dynamics. This opens the interesting possibility 
of doing quantum computation "straddling" a DFS, thus hav- 
ing a bigger code for computation, that still has much better 
protection than decoherence than the ones that are orthogo- 
nal to the DFS. Moreover, conditions for the existence of DFS 
are considered unstable. Our results show that an approximate 
notion of decoherence-free subspace is possible and useful for 
protecting quantum information. 

We also proved that this system thermalizes, but not always 



|2) = ^[-1,0,0,1], 



|3) = ^[0,-1,1,0], 



|4) = [r^,s^,s^,r^]. 
For notation simplicity we have introduced 



r± = Lll± 



and 



2V Vw^ + 4A2 



s± = ±A[4A2 + v{v ± a/w2+4A2)]-3. 



2. Matrix elements of the reduced density matrix 

We present here the exphcit calculation of the rates appear- 
ing in the Master equation |5] that we use to solve for the dy- 
namics of the reduced density matrix pg. As explained in the 
main text in the case of one single bath for both spins a (one 
dimensional) DFS exists. Because of the absence of bias terms 
in the Hamiltonian Hs this DFS happens to coincide with the 
eigenvector |3). 

Uncorrelated baths — Assuming no correlation between 
the two baths,according to I161I we have 



r+.,„ = jE("*l^'^l^)(^l^^l") 



X / dtexp{~iujint){Xi{t)Xi)i3, 

JO 



(12) 
(13) 



and 



K^kin = I E M^l\k){lK\n) (14) 



X / dtcxp{-iu;^kt){XiX,{t))p. (15) 

Jo 

Where si = a^ and s^ = r^, and X^ = J2a^aAbi,i + 
ba,i), we use the interaction representation Xi{t) — 
exp{—iHBt)Xi expiiHst). Here {■)p is the thermal equi- 
librium average over the bath's degrees of freedom. 
A direct calculation shows that 



{X^{t)X,) 



dujJi{uj) 



coth( ) cos(a;t) — i smiuit) 



The only quantity we need to compute is an integral of the 
form: 



^^M 



dt e~^"* 



duj' Ji{uj') 



coth(— — ) cos{uj't) =F i sin(w'i) 



For w > we obtain (the case w < is obtained by com- 
plex conjugation): 



It{u) = -Uu:) 



coth ^— T 1 
2 



iV 



duJ Ji{LJ ) 



bJ 



■ coth — =F 



LO 



V indicates the Cauchy principal value. The imaginary part 
can be absorbed as shift of the free oscillation frequencies of 
the system. For simplicity we disregard it. Because of the 
symmetries of Hs we only need to calculate W21 and W31, 
we obtain 



W. 



31 






coth( 



(3(^31 



and 



W2l = 



irA^k 
Vv^ + A? 



coth(^— ) - 1 



We report here the explicit solution of the diagonal ele- 
ments of the density matrix in the diagonal basis. We intro- 
duce Fi — W21 + W\i and Yi — Wj,\ + W^ia, using this 
notation we have 



Pii(t) 



l + e''"'^0 fl + e*^"-^' 



-t(ri+r2) 



+e' 



/3a;2i 



-l+e 



tri 



) (1 + e^"-) p22(0) + e^"- (-1 + e*^^) (e^^- (-1 + e'^^) + (l + e^"-) p33(0))] , 



P22(t) = [(1 + e""-) (1 + e^"-)] "' e-*(r^+r^) [(- 

'l + p/"^Mp22(0) + e 



gtr2+;3a.3i ^-y ^g/3u;2l\ _|_ gtri _|_ g/3("21+"3l) _|_gtri+/3 



+ e ^ + e "^^ 



li^-i\_ 



-1 + e^ 



^-l)pii(O) 



-e^"" (-1 + e*^^) (1 + e""-^!) ^22(0) + (e*^^ + e'^^^'^) (e^^^i (-1 + 6*^^) + (l + /'^'^) PssiO))] 

I 



The calculation of the off-diagonal elements does not main text. The explicit expressions of the 7's are directly re- 
present any difficulties and it obeys the relation given in the lated to W21 and W31 via Eq[8] 



9 



Single bath — In the case of one single bath we have 



states, i.e. W^i = W32 = W34 = as the direct calculation 
shows. In this case P33(i) = ^33(0). We have 



^t^kln = ^(mlcT;, +T^|A:)(/|cr^ +T^|n) 



and 



X / dteyi^{-iujint){X{t)X)p, 
Jo 



and 



Wi; 



Vv^ + A2 



coth( 



/3t^42> 



^mkin = l{M<^z + T,\k){l\cr^ + T^\n) 



X / dtexp{-iuJ„^kt){XX{t))lj. 
Jo 



W21 



2nA^k 

^y2 _|_ ^2 



coth(^-) - 1 



The two coupled equations for pu and P22 given by Eq. |6] 
As anticipated the state |3) is totally uncoupled to the other are solved by Laplace transform, we obtain 



Pii(A) 



A (A + W12 + W24 + W42) Pll(O) + W12 (Ap22(0) - W24 (/933(0) - 1)) 

A (1^12 (A + W24) + (A + T^2i) (A + W24 + W42)) 



P22(A) 



W2l{W24 + A(pii(0) + P22(0)) ~ Ty24Ap33(0)) + A(Ap22(0) ~ 1^24(^11(0) + Ap33(0) ~ 1)) 
A (W^12 (A + 1^24) + (A + W^2l) (A + VK24 + 1^42)) 



r 



The fourth component P44 is obtained from the normaliza- 
tion condition. To transform back to time domain the previous 
expressions does not comport any difficulties. Unfortunately 
the results can not be cast in a compact form so we do not 
report them here. Again the off-diagonal elements of the re- 
duced density matrix obey the expression given in the main 
text but in in this case all the rates 7's are expressed via 14^21 
and W42- We report the relevant ones here for convenience. 



713 = {W2i)/2 

723 = (M^12 + W42)/2 
743 = (M^24)/2 



(16) 
(17) 
(18) 
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